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INTRODUCTION 


A  number  of  techniques  are  currently  available  for 
seismic  discrimination  and  yield  estimation.  Most  of  these 
methods,  however  fail  to  fully  utilize  the  available 
azimuthal  information  and,  indeed,  many  methods  for 
magnitude  estimate  use  only  a  single  amplitude  measure. 
Shallow  earthquakes  will,  in  general,  produce  waveforms  and 
amplitude  patterns  that  vary  azimuthally  while  a  pure 
explosion  in  theory  produces  azimuthal ly  independent  wave¬ 
forms  and  amplitudes.  This  fact  is  not  utilized  in  most 
discrimination  schemes.  In  the  yield  estimation  context, 
azimuthal  variations  in  waveform  and  amplitude  may  contain 
information  on  tectonic  release  for  which  a  correction 
should  be,  but  is  usually  not,  made  in  obtaining  yield 
estimates.  In  this  document,  we  describe  a  computer 
program,  MOTIVE,  which  makes  use  of  this  azimuthal 
information,  as  well  as  utilizing  significant  portions  of 
the  waveform  information,  in  both  discrimination  and  yield 
estimation  contexts. 

MOTIVE  is  a  joint  body  wave  and  surface  wave  moment 
tensor  inversion  program.  Given  properly  windowed  body  wave 
and  surface  wave  seismograms  and  a  trial  depth,  MOTIVE 
determines  the  second  order  moment  tensor  that  provides  the 
best  fit,  in  a  least  squares  sense,  of  synthetics  to  the 
data.  By  using  a  number  of  trial  depths  and  selecting  the 
depth  that  produces  the  minimum  error,  depth,  as  well  as 
source  type  and  orientation,  may  be  determined. 

MOTIVE  incorporates  several  unique  features  in  addition 
to  the  joint  use  of  body  wave  and  surface  wave  data.  The 
user  can  specify  any  of  four  source  types;  an  unconstrained 
moment  tensor,  an  isotropic  source,  a  double  couple  source 
and  a  double  couple  plus  isotropic  source.  This  gives  the 
user  the  ability  to  explore  the  sensitivity  of  his  solution 
to  the  most  commonly  made  physical  assumptions  about  the 
source.  Moreover,  MOTIVE  can  include  state-of-the-art  body 
wave  and  surface  wave  path  and  receiver  corrections, 
determined  for  each  source-station  pair.  This  should  help 
eliminate  one  of  the  major  sources  of  error  and  bias  in 
moment  tensor  estimation,  particularly  where  short-period 
body  wave  seismograms  are  used. 

A  flexible  modular  approach  has  been  used  in  developing 
MOTIVE,  in  order  that  the  same  program,  through  choice  of 
user  options,  may  be  used  for  both  routine  processing  and 
research.  This  flexibility  also  insures  that  future 
improvements  can  be  readily  incorporated  into  the  existing 
program  structure. 
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ABSTRACT 

This  document  contains  the  functional  specifications 
for  a  program  to  invert  for  source  properties  using  a  moment 
tensor  source  description.  The  program  structure  is 
discussed;  individual  subroutines  are  named  and  their 
function  is  specified,  and  common  blocks  are  named  and  their 
variables  identified.  Though  not  a  final  working  program, 
these  functional  specifications  determine  the  orientation, 
flow  and  interaction  of  the  software  in  detail  which  will  be 
exceeded  only  by  the  code  itself. 
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PROGRAM  STRUCTURE 

The  moment  tensor  package  MOTIVE  may  be  conveniently 
divided  into  several  functional  modules,  each  containing 
several  subroutines.  Communication  between  modules  is 
through  common  blocks  and  data  files.  By  maintaining  the 
functional  independence  of  these  modules,  it  becomes  a 
relatively  simple  matter  to  incorporate  major  modifications 
in  the  overall  system,  both  in  a  research  mode  or,  if 
warranted  by  future  research,  in  a  production  mode.  In 
addition,  it  becomes  possible  to  maintain  several  options 
for  any  functional  portion  of  the  program,  where  option 
choice  is  specified  by  a  user  set  flag.  This  would  allow 
the  user  to  assess  the  effect,  for  instance,  of  a  different 
choice  of  error  function,  and  hence  different  partial 
derivatives,  on  the  final  solution,  without  a  major 
restructuring  of  the  program. 

A  flow  diagram  illustrating  all  major  program  modules 
is  shown  in  Figure  1.  The  first  module  encountered  is  the 
input  module.  This  module  consists  of  a  single  subroutine 
which  reads  station  independent  parameters,  such  as  source 
model,  the  number  of  body  and  surface  wave  stations  to  be 
employed,  which  program  options  are  to  be  used,  and  trial 
source  depths. 

The  body  wave  and  surface  wave  modules  each  ha^e  two 
basic  functions.  These  are  the  input  and  preprocessing  of 
the  appropriate  data  type  for  each  station  and  the 
generation  of  primitive  Green's  functions  and  partial 
derivatives  for  each  moment  tensor  element  and  station.  As 
the  body  and  surface  wave  inversions  employ  different  types 
of  data,  the  two  modules  must  have  somewhat  differing 
structures,  although  parallel  construction  has  been  used 
where  possible. 

The  data  used  by  the  surface  wave  portion  of  the 
inverse  are  complex  amplitudes  at  several  frequencies. 
Thus,  the  surface  wave  program  must  obtain  a  seismogram  from 
a  data  files  and  perform  a  Fourier  transform  (or  narrow  band 
filter)  to  obtain  a  spectral  estimate.  From  the  source  to 
station  distance  and  azimuth.  Green's  functions  at  the 
appropriate  frequencies  are  then  computed,  using  table 
lookup  for  a  standard  model.  These  Green's  functions  are 
then  corrected  for  the  specific  path,  where  such  information 
is  available.  These  corrections,  which  include  attenuation, 
dispersion,  multipathing  and  site  amplification  effects,  are 
to  be  provided  by  Systems,  Science  and  Software,  Inc.  in  the 
form  of  a  table  containing  corrections  for  each  source 
region-station  pair.  Data  and  Green's  function  spectral 
amplitudes  are  then  packed  into  data  and  partial  derivative 
arrays,  respectively,  for  use  by  the  inversion  module. 

Data  used  by  the  body  wave  portion  of  the  inverse  are 
in  the  form  of  time  series.  Thus,  for  each  station,  data 
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PROGRAM  MOTIVE 


Figure  1.  Flow  diagram  for  program  MOTIVE 
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must  be  obtained  from  a  data  file,  aligned  in  time, 
prefiltered,  windowed  and  normalized.  Green’s  functions  are 
computed  in  two  stages,  corresponding  to  depth  independent 
and  depth  dependent  portions.  The  depth  independent  portion 
consists  of  geometric  spreading,  residual  attenuation, 
receiver  function,  instrument  and  assumed  source  time 
function.  Of  these,  geometric  spreading  and  attenuation  are 
obtained  by  table  lookup  and  subroutine  calculation, 
respectively,  while  instrument  and  receiver  functions  are 
obtained  from  pre-existing  files.  These  files  are  accessed 
using  information  provided  on  input  for  each  station. 

The  depth  dependent  portion  of  the  Green's  functions 
are  computed  in  a  Haskell  matrix  program  for  each  moment 
ensor.  These  depth  dependent  portions  are  then  convolved 
with  the  depth  independent  preliminary  Green  functions.  The 
resulting  Green's  functions  are  time  aligned,  normalized, 
filtered  and  windowed  in  precisely  the  same  way  as  were  the 
data  for  that  station.  Data  and  Green's  fnctions  are  then 
packed  into  data  and  partal  derivative  matrices, 
respectively,  for  use  by  the  inversion  module. 

The  inversion  module  uses  the  partial  derivative  and 
data  arrays  constructed  in  the  body  and  surface  wave  modules 
to  construct  a  moment  tensor  that  minimize  the  difference 
between  model  and  data  in  a  weighted  least  squares  sense. 
Non-linear  constraints  are  iteratively  applied  to  a  linear, 
stabilized  least  squares  inversion,  so  that  optimal 
isotropic,  double  couple,  isotropic  plus  double  couple,  or 
unconstrained  moment  tensor  solutions  are  available  to  the 
user. 


Following  the  inversion  module  is  an  optional  time 
function  estimation  module,  which  will  be  discussed  below. 
The  optimal  moment  tensor  solution  for  that  depth,  together 
with  the  Green's  function  partial  derivative  arrays  and  the 
error  function  for  that  depth  are  passed  to  the  decision  and 
statistics  module.  Here  the  source  depth,  moment  tensor  and 
total  error  are  stored.  If  the  total  error  is  a  minimum 
compared  to  tht  determined  for  other  depths,  the  Green's 
functions  and  partial  derivatives  are  also  stored.  The 
program  then  loops  back  to  the  body  wave  module  in  order  to 
process  the  next  trial  source  depth. 

When  the  iterations  over  source  depth  are  completed, 
the  decision  module  outputs  a  table  of  depths,  moment  tensor 
solutions  and  total  errors .  For  the  depths  at  which  minimum 
error  occurs,  a  variance  and  data  importance  matrix  are  also 
produced.  For  comparison  purposes,  at  the  option  of  the 
user,  plots  of  the  solution  time  function  may  be  obtained, 
if  this  option  has  been  implemented,  as  well  as  plots  of 
data  and  synthetic  seismograms. 


CALLING  SEQUENCE 

MOTIVE 

FLAP 

MOD  IN 

SENSE 

SURFIN 

SEISIN 

UPNORM 

COOLB 

GBANDP 

PAW 

SEGRIN 

SWCSSS 

UPNORM 

GBANDP 

PAW 

PENSE 

RAYPRM 
POD  IN 

SEISIN 

SMOOTH 

UPNORM 

PAW 

PGR  IN 

Q2 

GEOM 

RRFRD 

RDINST 

HACK 

REMOD 

WRATH 

HASK2 

SMPHS 

INTRP 

GRAFT 

COOLB 

SMOOTH 

TSHFT 

UPNORM 

PAW 

MOTEST 

ASETUP 

MATACM 

PACK2 

WT2 

VECMAT 

MATINV 

MINV 

LSOLVE 

MSPLIT 

EIGEN 

PACK 

VC2ANG 
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Program  MOTIVE 

This  is  the  main  program  controlling  input  of 
data,  computation  of  Green's  functions,  inversion  for 
moment  tensor  components  and  output  of  results.  The 
program  is  executed  by  seven  modules,  three  of  which 
are  basically  input,  three  of  which  are  executed  in  a 
loop  on  depth,  and  the  last  one  of  which  is  a  summary 
and  output  module. 

CALLS:  (1)  FLAP  -  Reads  station  independent 

flags  and  parameters.  Surface 
wave  input  and  set-up.  Reads 
and  preprocesses  surface  wave 
data.  Defines  surface  wave 
Green's  functions.  Makes  both 
available  to  moment  tnesor 
inversion  module. 

(2)  SENSE  Surface  wave  input  and  set-up. 

Reads  and  preprocesses  surface 
wave  data.  Defines  surface 
wave  Green's  functions.  Makes 
both  available  to  moment 
tensor  inversion  moduel . 

(3)  PENSE  -  P-wave  input  and  set-up. 

Reads  and  preprocesses 
observed  P-wave  seismograms 
and  packs  data  into  data 
vector  with  surface  waves. 
Computes  preliminary  P-wave 
Green's  functions  lacking  only 
the  source  crustal  interaction 
which  must  be  done  in  loop  on 
depth. 

(4)  HACK  -  Thompson-Haskell  matrix 

package  to  compute  source 
crustal  response  as  a  function 
of  source  depth  and  moment 
tensor  component. 

(5)  GRAFT  -  Computes  final  P-wave 

Green ' s  functions  for  each 
moment  tensor  component  for 
the  current  depth  and  packs  in 
same  Green's  function  vector 
with  surface  wave  Green's 
functions . 
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(6)  MOTEST- 


(7)  DECIDE - 


Moment  tensor  estimation 
Performs  the  moment  tensor 
inversion  to  final  best 
isotropic  source,  best  double 
couple  source  or  best 
isotropic  plus  double  couple 
source,  as  requested,  for  the 
current  depth. 

Decision  and  statistics 
module.  Compares  results  of 
moment  tensor  inversion  at  the 
depths  tested;  determines 
overall  best  fit  and  computes 
supporting  statistics.  Does 
graphic  and  tabular  output  of 
final  results. 


Input: 

Output: 


MOTIVE  does  no  direct  input 
MOTIVE  does  no  direct  output 


Subroutine  FLAP 


Does  flag  and  parameter  input,  and  inputs  the 
source  crustal  model.  Flags  input  here  are  station 
independent. 

Called  by:  MOTIVE 


Calls:  MODIN 

Input:  none 

Reads :  NSTAS ,  NSTAP 

DTS ,  DTP 

NPS,  NPP 

TSMS ,  TSMP 

HMIN 

DH 

NH 

NITER 

ITYPE 

option  flags 


reads  crustal  model  file. 


numbers  of  stations  for 
input  of  surface  wave  and 
P-wave  data,  respectively 

sample  intervals 

desired  data  lengths 

smoothing  lengths 

minimum  source  depth 

increment  of  source  depth 

number  of  source  depths 

number  of  interations  in 
moment  tensor  inversion 
module 

determines  source  type  in 
the  moment  tensor 
inversion 

currently  undefined 


Output  flags  are  made  available  to  other  routines  through 
common  blocks. 
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Subroutine  MODIN 

Routine  reads  name  of  diskfile  containing  source 


crustal  model. 
Model  data  are 
in  COMMON/MOD/ 

opens  file,  reads  data,  and  closes  file, 
made  available  to  Haskell  matrix  routine 

Called  by  : 

FLAP 

Calls  : 

none 

Input  : 

none 

Reads  : 

NLAY 

number  of  layers  in 
crustal  model 

ALPHA 

layer  P-wave  velocity 

BETA 

layer  S-wave  velocity 

DENS 

layer  density 

THICK 

Layer  thickness 

Output:  Read  values  are  made  available  to  other  routines 

in  COMMON/MOD/. 


Subroutine  SENSE 


This  is  the  major  calling  subroutine  for  entry 
of  observed  surface  wave  data  and  generation  of 
corresponding  Green's  functions. 


Called  by  : 

MOTIVE 

Calls  : 

SURFIN 

Inputs,  manipulates  and  stores 
observed  surface  wave  data  in 
parameteric  form 

SGRIN 

Generates  surface  wave  Green's 

functions  and  stores  parametric 
representation  corresponding 
to  observed  data. 


Input  :  NSTAS  (in/FLAGS/)  -  number  of  stations  with 

surface  wave  data 

Reads  :  STID  -  4  column  station  indentifier 

DEL  -  epicentral  distance  in  degrees 

A2  -  azimuth 

QEFF  -  effective  Q  (optional) 

WT  -  estimate  of  station  or  record 

quality 


Output 


None 


Subroutine  SURFIN  (KS) 


This  is  the  calling  subroutine  for  entry  and 
manipulation  of  observed  surface  wave  data.  Input  data 
is  parameterized  by  its  complex  amplitude  at  selected 
frequencies  before  storage  in  a  data  vector. 


Called  by  : 

SENSE 

Calls  : 

SEISIN 

- 

reads  seismogram  from  diskfile 

UPNORM 

- 

normalizes  seismogram  to  unit 
power 

COOLB 

- 

utility  Fast  Fourier  Transform 
routine 

GBANDP 

- 

Gaussian  bandpass  filter 

PAW 

- 

Puts  parameterized  seismogram 
into  data  vector 

Input  : 

KS 

station  number 

Output 

FNORM  - 

normalizing  factor  for  current 
seismogram 

DV 

data 

vector  output  through  PAW 
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Subroutine  SEISIN  (SEIS,  NPP) 

This  routine  reads  window  flags,  file  name,  opens  disk 
file,  reads  seismogram,  closes  disk  file  and  puts  seismogram 
in  proper  form  for  later  manipulations. 

Called  by:  SURFIN,  SGRIN,  PODIN,  GRAFT 

calls:  OPENIN  -  open  a  disk  file  for  input 

RDDSK  -  system  routine  for  direct  disk  I/O 
CLOSE  -  closes  diskfile 
INTERP  -  linear  interpolator 


Input: 

NPP  - 

desired  number  of  points  for  P-wave  seismo¬ 
gram 

Output: 

SEIS  - 

seismogram 
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Subroutine  UPNORM  (X,  NP,  SNORM) 

Determines  total  power  in  array  X  and  normalizes  array 
X  to  unit  power.  If  SNORM  is  negative  upon  entry,  then  the 
routine  does  not  compute  power  but  uses  the  absolute  value 
of  SNORM  as  the  normalizing  factor. 

Called  by:  SURFIN,  SGRIN,  PODIN,  GRAFT 

Calls:  None 

Input:  X  -  array  to  be  normalized 

NP  -  length  of  X 

SNORM  -  power  in  array  X  upon  input 
Output:  X  -  normalized  to  unit  power 


Subroutine  COOLB  (NN,  XX,  SIGNI) 

Utility  subroutine  for  computation  of  Fast  Fourier 
Transforms . 

Called  by:  Various 

Calls:  None 

Input:  XX  -  complex  array  to  be  transformed 

NN  -  power  of  2  describing  length  of  XX 
SIGNI  -  Forward  transform  done  when  SIGNI  =  -1.0 
inverse  transform  done  when  SIGNI  =  1.0. 

Output:  XX  -  transformed  version 
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Subroutine  GBANDP  (X,  IF1,  IF2,  IF3,  IF4, ) 

Routine  to  perform  a  Gaussian  bandpass  operation 
on  the  spectrum  of  X. 

Called  by  :  SURFIN,  SGRIN 

Calls  :  None 

Input  :  X  -  Complex  frequency  domain  array 

:  IF1-IF4  -  frequency  designations  for 
limits  of  *•':  passband 

output  :  x  -  filtered  » 


< 
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Subroutine  PAW  (X,  NP,  IDG) 

This  routine  is  a  versatile  packing  program  which 
inserts  data  of  (potentially)  variable  lengths  into  a  single 
vector.  The  vector  type  is  data  vector  (IDG.LE.O)  or 
Green's  function  vector  (IDG.GT.O). 

Called  by:  SURFIN,  SGRIN,  PODIN,  GRAFT 

Calls:  None 

Input:  X  -  array  to  be  packed 

NP  -  length  of  X 

IDG  -  flag  determining  into  which  vector  X  will 
be  packed. 

Output:  ISTA,  JSTA  -  counter  of  how  many  X  arrays  have 

been  packed  in  vectors 

I COUNT,  J COUNT  -  counter  of  how  many  points  have 
been  packed  in  vectors 

DV  -  data  vector  of  length  ICOUNT 

containing  ISTA 
seismograms 

GV  -  Green’s  function  vector  of  length 

JCOUNT  containing  JSTA  Green's 
functions 
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Subroutine  SGRIN  (KS) 

This  is  the  calling  routine  for  calculating  surface 
wave  Green's  functions.  The  Green's  function  may  be 
visualized  as  the  impulse  response  of  the  total  travel 
path  from  source  to  receiver,  including  instrument.  The 
surface  wave  path  corrections  will  be  provided  by 
Systems,  Science  and  Software,  Inc.  The  Green's 
functions  are  parameterized  in  the  same  way  as  the  data 
and  are  stored  in  a  Green's  function  vector. 


Called  by  :  SENSE 

Calls  :  SWCSSS 

RDINST 

UPNORM 

COOLB 

GBANDP 

PAW 


Input  :  KS 

Output  :  GV 


Reads  surface  wave  path 
correction 

Reads  instrument  response 
and  returns  filter  for 
convolution 

Normalizes  Green's  function 
with  same  factor  as  used  in 
data 

Utility  Fast  Fourier  Transform 
routine 

Gaussian  bandpass  filter 

Packs  parameterized  Green's 
function  for  each  moment  tensor 
component  into  Green's  function 
vector 

station  number 

Green's  function  vector  output 
via  PAW 
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Subroutine  SWCSSS  (KS) 

Routine  reads  surface  wave  path  correction  function 
via  table  look-up  from  data  provided  by  Systems, 

Science  and  Software,  Inc.  The  path  correction  depends 
function  should  include  attenuation,  dispersion, 
multipathing  and  site  amplification  effects  and  should 
be  supplied  for  each  source  region-station  pair. 

Called  by  :  SGRIN 

Calls  :  undetermined 

KS  -  station  number 

SWGF  -  surface  wave  Green's  function 

in  the  frequency  domain 


Input 

Output 


Subroutine  PENSE 


This  is  the  major  calling  subroutine  for  entry  of 
observed  P-wave  data  and  generation  of  P-wave  Green's  func¬ 
tions  . 


Called  by. 

MOTIVE 

Calls : 

POD IN  - 

Inputs,  manipulates  and  stores 
observed  P  wave  data 

PGRIN  - 

generates  preliminary  P  wave 

Green's  functions  for  later 
convolution  with  Haskell  matrix 
response . 

Input:  NSTAP 

(in/FLAGS/) 

number  of  stations  with  P-wave 
data 

Reads:  STID 

_ 

4  column  station  indentifier 

DEL 

- 

epicentral  distance  in  degrees 

AZ 

- 

azimuth 

TSTAR 

- 

travel  time/effective  Q 

TAU2 

- 

absorption  band  parameter 

WT 

- 

estimate  of  station  or  record  quality 

Output:  Variables  read  in  are  passed  in  COMMON/STA/ 
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Subroutine  RAYPRM  (DEL,  P) 

Computes  ray  parameter  using  distance  and  a  table 
look-up  based  upon  the  Jeffries  -  Bullen  travel 
time  tables. 


Called  by  : 

PENSE 

Calls 

none 

Input  : 

DEL  - 

epicentral  distance  in  degrees 

Output  : 

P  - 

ray  parameter  in  COMMON/STA/ 
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Subroutine  PODIN  (KS) 

This  is  the  calling  subroutine  for  entry  and  manipula¬ 
tion  of  observed  P-wave  data  and  for  packing  of  that  data 
into  a  data  vector  for  use  in  the  inversion  module. 


Called  by:  SURFIN  PENSE 

Calls:  SEISIN  -  reads  seismogram  from  diskfile 

SMOOTH  -  performs  a  running  average  smooth¬ 
ing  of  the  seismogram.  (May  be 
replaced  by  a  bandpass  filter 
operation. ) 

UPNORM  -  normalizes  seismogram  to  unit 
power 

PACK  -  packs  seismogram  in  data  vector 
Input:  KS  -  station  number 


Output:  FNORM  -  normalizing  factor  for  current  seis¬ 

mogram.  Will  be  used  to  normalize 
Green's  function. 
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Subroutine  SMOOTH  (X,  NP,  DT,  TSM ) 

Performs  a  running  average  smoothing  of  array  X.  The 
number  of  points  in  the  running  average  operator  is  defined 
by  TSM/DT . 

Called  by:  PODIN,  GRAFT 

Calls:  None 

Input:  X  -  array  to  be  smoothed 

NP  -  number  of  points  in  X 
DT  -  sample  interval  for  X 
TSM  -  time  width  of  smoothing  operator  in  seconds 

Output:  X  -  smoothed  version 
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Subroutine  PGRIN  (KS) 

This  is  the  calling  subroutine  for  the  calculation  of 
preliminary  Green's  functions  for  P  waves.  The.  Green's 
function  may  be  visualized  as  the  impulse  response  of  the 
total  travel  path  from  source  to  receiver,  including  instru¬ 
ment.  This  routine  calculates  the  response  of  all  parts  of 
the  path  except  the  source  crustal  response  which  must  be 
done  in  an  iteration  on  source  depth. 


Called  by: 

PENSE 

Calls : 

Q2 

GEOM 

RDRRF 

RDINST 

computes  the  anelastic  attenuation 
computes  geometric  attenuation 
reads  relative  receiver  function 
and  prepares  for  convolution 
reads  instrument  response  and 
prepares  for  convolution 

Input: 

KS 

station  number 

Output: 

PWC 

P-wave  Green's  function.  Actually 
a  pre¬ 
liminary  or  primitive  Green's 
function 

in  that  neither  the  source  nor  the 
source 

crustal  response  are  included, 
(output  in  COMMON/PGREEN/) 
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Subroutine  Q2  (DT,  TSTAR,  TAU2 ,  QF) 

Generates  the  analytic  transform  of  a  frequency  de¬ 
pendent  Q  filter  as  defined  by  Minster  (  Geophys  J.  R. 
Astr.  Soc.,  52,  P.  503,  1978). 

Called  by:  PGRIN 

Calls:  None 

Input:  DT  -  sample  interval  in  time  domain 

TSTAR  -  travel  time/affective  Q 
TAU2  -  absorption  band  parameter.  To  get  frequency 
independent  Q,  set  TAU2=  .001 

Output:  QF  -  the  desired  Q  filter  in  the  frequency  domain. 


v 
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Subroutine  GEOM  (DEL,  GSPRD) 

Computes  geometric  spreading  factor  based  upon  surface- 
to-surface  travel  path.  This  gain  factor  is  based  upon  a 
parameterization  of  a  curve  defined  by  Langston  (Ph.D. 
thesis.  Cal.  Tech.,  1978). 

Called  by:  PGRIN 

Calls:  None 

Input:  DEL  -  epicentral  distance  in  degrees 

Output:  GSPRD  -  the  geometric  attenuation  amplitude  loss 

factor 
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Subroutine  RRFRD  (RF) 

The  routine  reads  a  file  name,  opens  the  disk  file, 
reads  a  relative  receiver  function  in  the  time  domain, 
closes  the  file  and  manipulates  the  receiver  function  for 
convolution . 

Called  by:  PGRIN 

Calls:  OPENIN  -  opens  diskfile  for  reading 

CLOSE  -  closes  diskfile 

COOLB  -  Fast  Fourier  Transform  routine 
SPLIN3  -  spline  fit  interpolator 

Input:  None 

Reads:  RRFID  -  name  of  diskfile  containing  desired 
receiver  function 

Output:  RF  -  the  desired  relative  receiver  function 

in  the  frequency  domain 


29 


Subroutine  RDINST  (FI) 

This  routine  reads  an  instrument  response  file  name, 
opens  the  file,  reads  the  instrument  response  in  the  time 
domain,  closes  the  file  and  prepares  the  instrument  response 
for  convolution. 

Called  by:  SGRIN,  PGRIN 

Calls:  OPENIN  -  opens  disk  file  for  reading 

CLOSE  -  closes  disk  file 
COOLB  -  Fast  Fourier  Transform  routine 

Input:  None 

Reads:  INSTID  -  name  of  diskfile  containing  desired 

instrument  response 

Output:  FI  -  the  desired  filter  for  instrument  response 

in  the  frequency  domain 
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Subroutine  HACK  (KH,  KS ) 

This  is  the  calling  routine  for  the  Haskell  matrix 
module.  Computes  Haskell  matrix  response  for  the  six 
independent  moment  tensor  components .  This  routine  is 
called  inside  of  nested  loops  on  depth  and  station. 


Called  by: 

MOTIVE 

Calls : 

REMOD  - 

redefines  crustal  model  to  include 
artificial  layer  at  current  source 
depth 

WRATH  - 

Wave  response  of  Thompson  Haskell 
matrix.  Does  actual  computation 
of  the  response  functions 

I nput : 

KH  - 

depth  counter 

KS  - 

station  number 

HMIN  - 

(from  COMMON/DEPTH/)  minimum  depth 

DH  - 

(from  COMMON/DEPTH/)  depth 
increment 

Output:  None 
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Subroutine  REMOD  (H) 


Redefines  crustal  model  to  include  an  artificial  layer 
at  current  source  depth. 


Called  by: 
Calls : 
Input : 


Output: 


HACK 

None 

H  - 

current  source 

depth 

NLAY  - 

(from  COMMON/MOD/)  number 

of  layer 

ALPHA  - 

II 

velocity 

layer 

P-wave 

BETA  - 

II 

velocity 

layer 

S-wave 

DENS  - 

II 

layer 

density 

THICK  - 

II 

layer 

thickness 

(All  in  COMMON/MOD2/ ) 


NL  -  new  number  of  layers  =  NLAY  +  1 
VP  -  layer  P-wave  velocity 

VS  -  layer  S-wave  velocity 

RHD  -  layer  density 

TH  -  layer  thickness 

LIS  -  layer  including  source 


I 


l 


C 
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Subroutine  WRATH  (KS) 

Routine  computes  the  source  crustal  P-wave  response 
using  a  Thompson-Haskell  matrix  method  (see,  e.g.,  Haskell, 
J .  Geophys .  Res . ,  67,  4751,  1962).  The  layer  matrices  are 
manipulated  to  compute  separate  responses  for  each  of  the 
six  independent  moment  tensor  components. 


Called  by: 

HACK 

Calls : 

HASK1  - 

SMPHS  - 
INTRP  - 

computes  layer  matrices  at  a 
given  frequency 

unwinds  phase  of  response  function 
linear  interpolator  used  to 
expand  response  defined  at  limited 
set  of  frequencies  to  full 
frequency  range  for  convolution 

Input: 

KS  - 
/MOD2/  - 

P  - 

station  number 

see  definitions  in  subroutine 

REMOD 

(from  COMMON/STA/)ray  parameter 

Output:  HMR  (in  COMMON/HAS/) 

-  contains  the  response  functions 
for  the  six  tensor  components. 
This  doubly  dimensioned  array 
could  be  visualized  as  including 
HMRXX,  HMRYY,  HMRZZ ,  HMRXY, 

HMRXZ  and  HMRYZ 
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Subroutine  GRAFT  (KS) 


This  is  the  calling  routine  for  computation,  manipulation 
and  storage  of  final  P-wave  Green's  function  corrected  for 
Haskell  matrix  response. 


Called  by: 
Calls : 


Input: 


Output: 


MOTIVE 

Fast  Fourier  Transform  routine 
running  average  smoothing  operator 
apply  a  time  shift  to  give 
proper  alignment  to  Green's 
function 

normalize  to  unit  power 
put  current  six  components  of 
Green's  tensor  response  in 
Green's  function  vector. 

KS  -  station  number 

PWC  -  (in  COMMON/PGREEN/ )  P  wave 

correction  factor 

HMR  -  (in  COMMON/HAS/)  Six  Haskell  matrix 
responses,  one  for  each  moment 
tensor  component 

GV  -  Green's  function  vector  output  via 
PACK 


COOLB 
SMOOTH  - 
TSHFT  - 


UPNORM  - 
PAW  - 


> 
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Subroutine  TSHFT  (X,  DT,  SEC) 

Circular  shifts  X  array  by  time  'SEC'. 

Called  by:  GRAFT 

Calls:  None 

Input:  X  -  array  to  be  shifted 

DT  -  sample  interval 

SEC  -  amount  of  shift.  SEC  greater  than  zero 
causes  a  right  shift. 

Output:  X  -  shifted  version 
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Subroutine  MOTEST  (MODE,  MXLOOP,  RESTES,  MT) 

Principal  control  subroutine  for  moment  tensor  inversion 


Called  by: 

MOTIVE 

Calls : 

ASETUP  - 

Sets  up  angle  constants 

MATACM  - 

Accumulates  design  matrix  from 
each  data  point  matrix 

VECMAT  - 

Converts  vector  of  n(n-l)/2 
elements  to  n  x  n  matrix 

MATINV  - 

Matrix  inverter 

LSOLVE  - 

Matrix  multiplier 

MSPLIT  - 

Converts  best  fit  solution  matrix 
into  principal  values  and  returns 
appropriate  moment 

YPERTB  - 

Calculates  new  residuals 

RSTEST  - 

Tests  residual  and  number  of 
iterations  against  preset  values 

DATACM  - 

Finds  the  new  residual  vector 

CNSTR  - 

Determines  array  which  characterizes 
constraints  for  moment  tensor  matrix 

CONCOR  - 

Determines  the  matrix  which  relates 
the  constrained  and  free  solutions 

VARNC  - 

Calculates  95%  confidence  intervals 

Input: 

MODE 

Solution  type 

1.  Free  (no  constraints) 

2 .  Double  couple 

3 .  Double  couple  and  isotropic 

4.  Isotropic 

MXLOOP  - 

Maximum  number  of  iterations 

RESTES  - 

Residual  test  level 

Output: 

MT  - 

Best  estimate  of  moment  tensor 

S* '  r 
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Subroutine  ASETUP 


Sets 

Called  by: 
Calls : 
Input: 
Output : 


» 


I 


c 


up  angle  constants  (radians  to  degrees,  etc. 
MOTEST 
None 
None 

Passes  arguments  in  labeled  common  - 
COMMON/CONST/P I,  R2DEG,  DEG2R ,  SQRT2 
PI  :  n 

R2DEG:  conversion  factor  for  radians 

to  degrees 

DEG2R :  conversion  factor  for  degrees 

to  radians 

SQRT2 :  square  root  of  2 
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Subroutine  MATACM  (M,  FR,  NQ,  YR,  XM,  YVEC ) 

Accumulates  design  matrix  XM  and  data  vector  YVEC  from 
each  data  point. 

Called  by:  MOTEST 

Calls:  1)  PACK2  -  Converts  from  natural  elements  of 

the  moment  tensor  to  double  couple  and 
isotropic  form 

2 )  WT2  -  determines  the  weighting  factors 

INPUT:  M  :  number  of  data  .. 

FR:  vector  of  partials  £or  in  datum 

NQ:  quality  factor  of  nr  datum 
YR:  rrt  datum 

T  th 

Output:  XM:  ( FR )  V  ( FR ) ;  F  is  the  matrix  of  m  partial,  V 

is  the  variance  matrix 
YVEC:  (FR)T  V  (YR) 


1 


I 
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Subroutine  VECMAT  (N,  FVEC,  FMAT ) 

Converts  FVEC,  a  vector  of  n(n-l)/2  elements  into  an 
n  x  n  matrix. 

Called  by:  MOTEST 

Calls:  None 

Input:  N  -  dimension 

FVEC  -  elements  of  (FR)1  V  (FR)  (Matrix  FMAT 

from  MATACM ) 

Output:  FMAT  -  matrix  representation  of  FVEC 


•  r* 
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Subroutine  MATINV  (FMAT,  NL) 

Inverts  FMAT  using  IBM  SSP  Gauss  Jordan  routines 


Called  by. 

MOTE ST 

Calls : 

MINV  - 

Matrix  inversion  routine  (SSP) 
Gauss  Jordan  method 

using 

Input: 

FMAT: 

the  NL  x  NL  design  matrix 

NL: 

dimension  of  FMAT 

Output : 

FMAT: 

inverted  form  of  input  design 

matrix 
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Subroutine  LSOLVE  (FMAT,  YVEC,  MSTART,  NL) 

Multiplies  inverted  design  matrix  FMAT  and  the  residual 
vector  YVEC. 

Called  by:  MOTEST 

Calls:  None 

Input:  FMAT:  NL  x  NL  design  matrix 

YVEC:  NL  x  1  residual  vector 

NL:  dimension 

Output:  MSTART:  NL  x  1  vector,  best  fit  estimate  of  the 

elements  of  the  moment  tensor 
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Subroutine  MSPLIT  (MSTART,  MOUNT,  PRIN,  I TYPE ) 

Converts  MOTEST  (best  estimate  of  the  elements  of  the 
moment  tensor)  into  its  principal  values  and  returns  to 
moment  tensor  conforming  to  the  mode  or  type  of  solution 
specified  in  MOTEST. 


Called  by: 

MOTEST 

Calls: 

EIGEN  - 

IBM  SSP  routine;  rotates  matrix  to 
determine  principal  components 

PACK 

Converts  elements  of  MSTART  to  or 
from  packed  mode  (double  couple  and 
isotropic) . 

VC2ANG  - 

Converts  vector  to  coordinate 
angles 

Input : 

MSTART  - 

Matrix  containing  the  6  elements 
of  the  best  fit  solution. 

I TYPE  - 

mode  of  solution 

Output : 

PRIN 

Eigenvectors,  eigenvalues  of  fJ 
[PRIN  (1)  =  1st  eigenvalue,  PRIN  (2 
and  PRIN  (3)  =  0,  4>  of  eigen¬ 
vector;  similarly  for  PRIN  (n), 
n  =  4,  5,  6,  7,  8,  9;  PRIN  (10)  = 
trace  (isotropic  part).] 

MOMNT  - 

Appropriate  typed  moment  tensor 
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Subroutine  RSTEST  (YVS,  LOOP,  NFLAG ) 

Tests  residual  and  number  of  iterations  against 
specified  values. 

Called  by:  MOTEST 

Calls:  None 


Input:  YVS  -  residual 

LOOP  -  number  of  the  iteration 

Labeled  COMMON  -  COMMOD/INITVL/MXLOOP , 

RESTES 

(max.  #  loops,  target 
residual  value  from 
MOTEST ) 

flag  to  indicate  whether 
to  exit  from 

iteration  scheme  (NFLAG  = 
2,  exit) 


Output: 


NFLAG 
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Subroutine  YPERTB  (MOMNT,  YVS ,  LOOP) 

Calculates  new  residuals  Y2  (Y2  =  Y  .  -  F-M)  to  be 

stored  in  labeled  common. 


Called  by:  MOSTART 

Calls:  WT2  -  determines  the  weighting  factor  V  for 

the  nth  datum. 


Input : 

MOMNT  - 

current  best  solution  for  the 
moment  tensor 

LOOP  - 

number  of  iterations 

Output: 

YVS 

residual 

Y2 

residual  matrix;  contained  in 
labeled  common:  COMMON/RDATA/ 
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Subroutine  RSTEST  (YVS,  LOOP,  NFLAG) 

Tests  residual  and  number  of  iterations  against 
specified  values. 

Called  by:  MOTEST 

Calls:  None 

Input:  YVS  -  residual 

LOOP  -  number  of  the  iteration 

Labeled  COMMON  -  COMMOD/INITVL/MXLOOP.RESTES 

(max.  #  loops,  target 
residual  value  from 
MOTEST ) 

Output:  NFLAG  -  flag  to  indicate  whether 

to  exit  from 

iteration  scheme  (NFLAG  = 

2,  exit) 


• 


-TB**  J 
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Subroutine  DATACM  (YVEC) 

Finds  the  new  residual  vector  [FTV  (Y0jjS  “  F*M)  =  YVEC] 

Called  by:  MOTEST 

Calls:  None 

Input:  (From  labeled  COMMON: 

COMMON/COEFF/N ,  F  (6,100),  F2  (6,100) 

COMMON /RDAT A/Y 1(400),  Y2(400) 
COMMON/QUAL/NQUA( 400 ) ,  V(400),  V2(400) 

V  -  variance  vector 

Y  -  observed  data  vector 

F  -  partial  derivative  matrix 

Output:  YVEC  -  New  residual  vector 
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Subroutine  CONSTR  ( ITYPE ,  NL,  D,  MOMNT) 


Determines  the  NL  x  6  matrix  which  characterizes  the 
constraints  for  the  moment  tensor 


Called  by: 

MOTEST 

Calls : 

None 

Input: 

MOMNT  - 

Current  best  estimate 
moment  tensor 

of  the 

ITYPE  - 

moment  tensor  type  (= 
from  MOTEST) 

MODE 

Output: 

NL  - 

#  of  lines  in  D 

D  - 

constraint  matrix 
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Subroutine  CONCOR  ( FMAT ,  D,  NL,  ,  B) 

Determines  the  matrix  that  will  relate  the  constrained 
and  free  solutions. 


Called  by: 

MOTEST 

Calls : 

MAT I NV  - 
MAP ROD  - 

matrix  inversion  routine 
matrix  multiplication  routine 

Input : 

FMAT 

D 

NL 

design  matrix  (the  free  solution) 
constraint  matrix 
#  of  lines  in  D 

Output: 

B 

NL  x  NL  matrix  relating  free  and  con¬ 
strained  solutions. 
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Subroutine  VARNC  ( FMAT,  YVS,  B,  PRIN,  DPRIN ) 

Calculates  the  95%  confidence  intervals  of  the  principal 
components . 


Called  by: 

MOTEST 

Calls : 

HEXT  - 

Rearranges  inputs  for  VARNC 

Input: 

FMAT  - 

design  matrix 

YVS  - 

Design  matrix 

B  - 

matrix  relating  constrained  and 
free  solutions 

PRIN  - 

Eigenvalues  and  eigenvectors  of  the 
moment  tensor 

Output: 

DPRIN  - 

95%  confidence  intervals  of  PRIN 
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Subroutine  PACK2  (FIN,  FOUT) 

Converts  from  the  natural  elements  of  the  moment  tensor 
to  double  couple  and  isotropic  form.  This  solely  for  the 
partials  of  F  where  A  =  FIN  M  =  FOUT  M' 

A  =  amplitude 
FIN  =  original  factors 
M  =  natural  elements 
FOUT  =  new  factors 
M'  =  packed  elements 


Called  by: 

MATACM 

Calls : 

None 

Input: 

FIN  - 

natural 

elements  of  M 

Output: 

FOUT  - 

double 

couple  and  isotropic  factors 

of  M 


Subroutine  WT2  (M,  YR,  LOOP) 

Determines  the  weighting  factor  for  the  mth  datum. 
Arguments  are  passed  via  COMMON. 


Called  by: 

MATACM 

YPERTB 

Calls : 

None 

Input : 

M  - 
YR  - 

data  number 

Mth  residual 

Output : 

V  - 

V2  - 

variance  function  of 
variance  function  of 
and  quality  factor 

quality  factor 
the  residual 

(Output  via  COMMON/QUAL/ ) 
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Subroutine  PACK  (MOMIN,  MOMOUT,  JMODE) 

Converts  elements  of  moment  tensor  to  or  from  packed 

form . 

Called  by:  MSPLIT 

Calls:  None 

Input:  MOM IN  -  input  moment  tensor 

JMODE  -  flag  indicating  whether  to  convert 
MOM IN  to  (JMODE  -  1 )  or  from 
(JMODE  =  2)  packed  form. 

Output:  MOMOUT  -  Packed  or  unpacked  form  of  MOM  IN 


Subroutine  VC 2 AN G  (R,  THETA,  PHI) 
Converts  vector  to  angles 


0  increases  counter¬ 
clockwise  from  vertical 
down 


Called  by:  MSPLIT 

Calls:  None 

R  -  input  vector 


Input: 

Output: 


THETA 


0 


I 


i 

* 
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Subroutine  HEXT  (A,  EIG,  EIGV,  SIGMA,  C) 

Arranges  inputs  for  VARNC 
Called  by:  VARNC 

Calls:  None 

input:  A  -  matrix  to  be  analyzed 

EIG,  EIGV  -  principal  components 

v  -  covariance  matrix  of  A 
SIGMA  -  standard  error  of  residuals 
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Subroutine  MAPROD  (A,  B,  I,  J,  K,  C) 

Matrix  multiplier;  calculates  C  =  A  B  where  A  is  a  i  x  j 
matrix;  B  is  a  j  x  k  matrix;  C  is  a  i  x  i  matrix,  i,  j, 
and  k  5  6 

Called  by:  CONCOR 

Calls:  None 

Input:  matrices  A  and  B 

matrix  dimensions  I,  J,  K 

Output:  matrix  C 


1 


< 
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Subroutine  DECIDE  (H) 

Major  subroutine  in  decision  and  statistics 
module.  Saves  and  ultimately  prints  moment  tensor 
elements  and  total  error  for  each  depth.  Computes  data 
importance  matrix  for  depth  at  which  minimum  error 
occurs .  Prints  variance  matrix  for  that  depth. 


Called  by  : 

MOTIVE 

Calls 

OUTPUT 

- 

Does  plot  output 

Input  : 

H 

- 

depth 

GV 

- 

(in  COMMON/GREVEC/ )  Green’s 
function  vector 

DV 

- 

(in  COMMON/DATVEC/)  data 
vector 

Output 

GRSV 

- 

best  fit  Green's  function 
array 

SVAR 

- 

variance  matrix 

SAVSTF 

Depth,  moment  tensor  solution, 
amplitude  of  double  couple 
component  and  total  error  for 
each  depth. 

•>  ' 
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Subroutine  OUTPUT  (X,  Y,  PLTID) 

Plots  data  (X)  and  best  fit  synthetic  (Y). 
Uses  system  plotting  routines. 

Called  from  :  DECIDE 

Calls  :  system  dependent  plot  routines 

Input  :  X  -  data  seismogram 

Y  -  synthetic  seismogram 

PLTID  -  plot  indentifier 

Output  :  graphic  display  of  results 
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List  of  COMMON  blocks. 

Detailed  explanation  of  variable  names  are  given 
on  the  following  pages.  The  arrays  must  be  dimensioned 
to  be  equal  to  or  larger  than  minimum  sizes  determined 
by: 

NS  :  number  of  stations.  Counts  surface 

waves  and  body  waves  as  separate  stations 

NSS  :  number  of  stations  contributing  surface 

wave  data 

NSP  :  number  of  stations  contributing  body 

wave  data 

NH  :  number  of  depths 

NC  :  number  of  moment  tensor  components 

( always  =  6 ) 

NT  :  number  of  time  samples  (NPS  or  NPP) 

NF  :  number  of  frequency  samples  (always  512 

when  stored) 

NL  :  number  of  layers  in  crustal  model 

If  no  dimension  specifications  are  provided,  then 
variable  is  not  an  array. 


I 


I 


t 


I 


i 
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COMMON/FLAGS/NSTAS ,  DTS,  NPS,  TSMS,  NSTAP ,  DTP,  NPP,  TSMP 
NSTAS ,  NSTAP 


TSMS,  TSMP 
NPS,  NPP 

DTS ,  DTP 
Required  by: 


numbers  of  stations  with 
input  data  for  surface  waves 
and  P  waves,  respectively 

time  widths  for  running 
average  smoothing  operators 

desired  numbers  of  points  in 
in  input  data  and  computed 
Green's  functions 

sample  intervals 

MOTIVE 

FLAP 

SENSE 

PENSE 

POD  IN 

PGR  IN 

HACK 

WRATH 

GRAFT 

MOTEST 
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COMMON/STA/STID  (NS),  DEL  (NS),  AZ  (NS).  QEFF  (NS),  TSTAR  (NS), 
TAU2  (S),  WT  (NS),  FNORM  (NS),  P  (NS) 


STID 

- 

four  column  station  identifier 

DEL 

- 

epicentral  distance  in  degrees 

AZ 

- 

Azimuth  in  degrees 

QEFF 

- 

effective  Q 

TSTAR 

- 

travel  time/effective  Q 

TAU2 

- 

absorption  band  parameter 

WT 

- 

estimated  station  variance  (quality) 

FNORM 

- 

factor  required  to  normalize 
seismogram  to  unit  power 

P 

- 

ray  parameter 

Required  by:  SENSE 

PENSE 
POD  IN 
PGR  IN 
HACK 
WRATH 
DECIDE 
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COMMON/DATVEC/ISTA,  I COUNT,  DV  (NSS  *  NPS  +  NSP  *  NPP ) 


ISTA 


I COUNT 
DV 


Required  by: 


station  counter  giving  number  of 
stations  with  data  in  the  data 
vector 

point  counter 

data  vector  of  length  I COUNT  into 
which  the  data  from  ISTA  stations 
have  been  packed 

PACK 

MOTE ST 


COMMON/GREVEC/ ISTA, 

JCOUNT,  GV(NSS  *  NPS  +  NSP  *  NPP) 

JUSTA 

“ 

station  counter  giving  number 
of  stations  with  data  in  the 
Green’s  function  vector 

JCOUNT 

- 

point  counter 

GV 

Green’s  function  vector  of  length 
JCOUNT  into  which  the  Green's 
functions  for  JSTA  stations  has 
been  packed 

Required  by: 

PACK 

MOTEST 

C 


C 
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COMMON/DEPTH/KMIN,  DH,  NH 

HMIN  -  minimum  source  depth 

DH  -  increment  of  source  depth 

NH  -  number  of  source  depths  to  be 

tested 


Required  by: 


FLAP 

MOTIVE 

HACK 


COMMON/INV/NITER,  I TYPE 


NITER 

I  TYPE 


Required  by: 


number  of  iterations  to  be  executed 
by  inversion  module 

specifies  type  of  source  to  be 
sought  in  inversion 

FLAP 

MOTEST 


COMMON/HAS /HMR ( NF ,  NC ) 

HMR  -  six  components  of  Haskell  matrix 

response,  corresponding  to  the  six 
moment  tensor  components  (XX,  YY, 
ZZ,  XY,  XZ,  YZ) 

Required  by:  WRATH 


.  -r*  -’/T3" rj 


GRAFT 


COMMON/GREEN/PWGF  (NF,  NC,  NSP),  SVG F  (NF,  NC,  NSS) 


PWGF 

SVGF 

Required  by: 


COMMON/PGREEN/PWC 

PWC 

Required  by: 


P-wave  Green's  Function 
Surface  wave  Green's  Function 
SENSE 
GRAFT 


(NC,  NSP) 

P-wave  path  correction  factor 
including  all  factors  but  the 
Haskell  matrix  response 

PGRIN 


GRAFT 


COMMON/MOD/NLA Y , 

ALPHA 

(NL),  BETA  (NL),  DENS  (NL),  THICK  (NL) 

NLAY 

- 

number  of  layers  in  crustal  model 

ALPHA 

- 

layer  P-wave  velocity 

BETA 

- 

layer  S-wave  velocity 

DENS 

- 

layer  density 

THICK 

- 

layer  thickness 

Required  by: 

MOD  IN 

REMOD 

COMMON /MOD2 /NLAY I , 

VP 

(NL),  VS 

(NL),  RHO  (NL),  TH 

NLAY  I 

- 

new  number  of  layers 

VP 

- 

layer 

P-wave  velocity 

US 

- 

layer 

S-wave  velocity 

RHO 

- 

layer 

density 

TH 

- 

layer 

thickness 

LIS 

- 

layer 

including  source 

Required  by: 

REMOD 

WRATH 

COMMON/PAS S/DEPTH,  TENS  (6),  AIS,  ADC,  PDC  (3),  VAR  (36),  ER 


DEPTH  -  Depth 

TENS  -  Moment  tensor  elements 

AIS  -  Amp i tude  of  isotropic  component 

ADC  -  Amplitude  of  double  couple 

component 

DCD  -  Double  couple  orientation 

VAR  -  Variance  matrix 

ER  -  Total  error 


Required  by: 


MOTEST 

DECIDE 


COMMON/SGR/GRSV  (512,  50,  6),  SVAR  (6,6),  SAVSTF  (13,30) 

GRSV  -  Greens  function  array  of  best  fit 

solution 

SVAR  -  Variance  matrix 

SAVSTF  -  Depth,  moment  tensor  solution, 

amplitude  of  isotropic  component, 
amplitude  of  double  couple 
component,  double  couple  orienta¬ 
tion  and  total  error  for  each 
depth 


Required  by: 


DECIDE 


